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Abstract —A dynamic model of the self-vibration cycle in an engine with opposed cylinders and two 
pistons located on a common guide bar is studied. To each cylinder containing a working liquid, a pair 
of hydraulic accumulators is attached the upper part of which contains a gas heated in one hydraulic 
accumulator and cooled in the other. Each hydraulic accumulator contains a membrane separating the 
gas and the liquid. In this way, a pressure difference is created in the cylinders that acts on each of the 
pistons and moves them along with the guiding bar. 
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At the current stage of progress in the machine-building industry, development of cyclic engines with 
external heating (Stirling-type engines) is becoming increasingly important. The operation principle and a 
description of the thermodynamic characteristics of the Stirling cycle are presented in detail in [ 1-3]. How¬ 
ever, in these and other studies by Russian and foreign researchers, analysis is performed on the basis of the 
thermodynamic cycle alone, i.e., without taking into account the forces of inertia and resistance that act on 
the moving parts of the engine mechanism. Neglect of these factors results in a significant decrease in the 
accuracy of calculations involved in the analysis of motion and sometimes even to qualitative distortions of 
the ideal cycle. In [4], based on the combined analysis of the thermodynamic process and the methods of 
machine dynamics applied as an example to the Stirling free-piston engine, it was shown that the isochoric 
parts of the ideal cycle [1-3] may not be implemented. In [5], a similar analysis was made, taking into 
account the forward-and-backward motion of a piston and rotation of a working shaft driven by the piston 
via a crank-and-rod mechanism, for the Stirling engine with a crosshead mechanism. 

Let us consider a model dynamic of an engine with opposed cylinders [6] (Fig. 1), where 1 is the first 
hydraulic accumulator (heating), 2 is the regenerator of the first hydraulic accumulator, 3 is the second 
hydraulic accumulator (cooling), 4 is the third hydraulic accumulator (cooling), 5 is the regenerator of the 
second hydraulic accumulator, 6 is the fourth hydraulic accumulator (heating), 7 is the right hydraulic cyl¬ 
inder, 8 is the right piston, 9 is the driving mechanism, 10 is the guiding bar, 11 is the left piston, and 12 is 
the left hydraulic cylinder. 

The engine’s working chambers consist of the working chambers of cylinders 7 and 12 with guiding bars 
connected in opposite directions. The double-acting cylinder is characterized by different efficient areas of 
the piston and piston-rod (left and right) ends. On the piping of piston ends, heaters (heat exchangers) are 
installed and the piping of the piston-rod ends bears coolers (cooling radiators). In the middle part of the 
piping connecting the ends, regenerators 2 and 5 are installed. To create an excessive pressure, gas is 
pumped into the working chambers through filling valves. The peri¬ 
odic motion of the pistons is guided by drive mechanism 9. The effi¬ 
cient areas of piston and piston-rod ends are pairwise equal. Owing 
to this, it is possible to maintain overpressure when the working 
chambers of an engine are filled. During the filling operation, the 
pressure values in the piston and piston-rod ends are the same, since 
the chambers are connected with piping and the hydraulic resistance 
of the heater, cooler, regenerator, and connecting piping is insignifi¬ 
cant, so that the excessive pressure during filling does not create any 
unbalanced forces. The unbalanced forces are created during the 
engine operation owing to the difference of pressures in cylinders 7 
and 12 operating in opposite phases. Fig- 
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The opposed location of the cylinders makes it possible to implement a positive feedback between the 
cylinders' working processes. This circumstance ensures a unique pattern of the engine’s periodic processes 
without introducing any adjustment devices [7] and enables obtainment of a significant initial pressure in 
the working chambers. When the flywheel of the driving mechanism rotates clockwise, the volume of the 
piston chamber in cylinder 7 decreases and that in cylinder 12 increases. This change causes motion of the 
working medium along the circulation contour; the major part of the gas is heated by the heater and the total 
pressure in cylinder 7 increases and in cylinder 12 decreases. This process creates an increasing force acting 
in the direction from the right-hand side of the wheel toward its left-hand side (positive feedback between 
the working processes in the chambers) to the pistons' dead center. The flywheel mechanically passes the 
dead center and causes motion in the opposite direction. One can see that continuous motion of the wheel 
is maintained in absence of any valves. The working chambers in this design do not need to be air tightened, 
since the pressure values in the left-hand and right-hand chambers of the hydraulic cylinders are pairwise 
equal and the pistons are installed to maintain circulation of the working medium. In this case, only rods of 
hydraulic cylinders have to be tightened. 

Here, we present a rigorous and a simplified model of the dynamics of the analyzed engine model. 

The rigorous model enables calculation of steady-state and transient regimes by numerical methods 
(using the MathCAD software package). The approximate model enables derivation of closed analytical 
solutions for the steady-state regimes of the engine’s operations. 

When developing the rigorous and approximate models, the following assumptions are made: deforma¬ 
tion of the membranes in the hydraulic accumulators when the liquid is displaced with the gas or vice versa 
is approximated by their displacement upwards or downwards, respectively; the working liquid is ideal and 
uncompressible; and pressures p h p 2 , p 3 , and p 4 developed in the upper (gas filled) parts of the accumulator 
volumes are passed directly to the working surfaces of the pistons. 

By virtue of the first assumption, the membrane surfaces are assumed to be flat. 

Let us introduce the following notations: V\ , V 2 , V 3 , and V 4 arc the total volumes of the liquid con¬ 
tained in the left-hand and right-hand chambers of the hydraulic cylinders and the lower chambers of the 
hydraulic accumulators; V 1A , V 2/ \, V 3A , and V 4A are the current values of the volume of the fluid contained in 
the lower chambers of the hydraulic accumulators; V lc , V 2c , V 3 C , and V 4t are the current values of the volume 
of the fluid contained in the left- and right-hand chambers of the hydraulic cylinders; S A is the area of the 
hydraulic accumulator’s membrane; S c is the area of the hydraulic cylinder’s piston; h is the thickness of the 
hydraulic cylinder’s piston; / is the height of the chambers of hydraulic accumulators; H is the length of the 
chambers of hydraulic accumulators; x is the current value of the piston displacement measured from the 
left end of each hydraulic cylinder; m ps is the mass of each piston; m ld is the rod mass; m t is the total mass 
of the engine’s moving parts; y b y 2 , y 3 , and y 4 are the current values of the liquid level in each hydraulic 
accumulator; r 1 and r 2 are the regenerator parameters characterizing the rate with which the gas mass 
changes as a function of the pressure difference; 7j and T 2 are the temperature of heating and cooling of 
hydraulic accumulators, respectively; m u m 2 , m 3 , and m 4 are the current values of the gas mass in each of 
the four hydraulic accumulators; m d is the invariable mass of the gas contained in the “dead volume” (heater, 
regenerator, and cooler); and m 01 and m 02 are the total masses of the gas contained in the first and second 
and in the third and fourth hydraulic accumulators. 

The balance of volumes for the first hydraulic accumulator and the left-hand side of the first hydraulic 
cylinder has the form 


Via + ^ic — + x $c — V\- 


( 1 ) 


From (1), it follows that 


= (F!-x5 c )/ 5 a . 


( 2 ) 


In accordance with (2), the volume of the gas contained in the upper chamber of the first hydraulic accu¬ 
mulator may be found as 


El = ( l-y\)S A = IS A - Vi + xS c . 


(3) 
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Since t), > 0, formula (3) yields the boundary of the leftmost position of the piston, x>x ]m = (Vi - IS A )/S C . 
The Mendeleev-Klaiperon equation for volume v t has the form 


PlVi 


m x 

M 


RJu 


( 4 ) 


where is the specific gas constant and M is the molar mass of the working medium. 

Formulas (3) and (4) determine the current values of the pressure p , (x, m ,) acting from the left on the 
piston in the first cylinder, 


P\{x, m i) 


m l RT l 

IS A — V i + xS c 


( 5 ) 


The distance between the piston and the right-hand wall of a hydraulic cylinder is x 2 = H-h-x. Respec¬ 
tively, the equation of the balance of volumes (1) for the right-hand side of the first hydraulic cylinder may 
be written as y 2 S A + (H-h-x)S c = V 2 • We obtain, then, 


3T = 


V 2 - (H-h -x)S c 

si 


( 6 ) 


Taking into account (6), we determine the gas volume in the upper chamber of a hydraulic accumulator, 

v 2 = lS A -V 2 + (H-h)S c + S c x. (7) 


Since v 2 > 0, formula (7) yields the rightmost position of a piston x < x rm = \ V 2 -IS A — (H- h)S c /S c ], 
The Mendeleev-Klaiperon equation for volume v 2 has the form 


m, 


PiV 2 - M R n T 2- 


( 8 ) 


From (8), we determine p 2 . 


Pi(x, m 2 ) 


_ m 2 RT 2 _ 

lS A -V 2 + (H-h)S c -S c x 


( 9 ) 


where R = R^/M is the universal gas constant. 

The expressions for the pressure values in the left- and right-hand chambers of the right-hand hydraulic 
cylinder may be derived similarly: 


p 2 = (x,m 3 ) 


m 3 RT 2 

1S A — V 2 + xS c 


Pa{x, m 4 ) 


m 4 RT l 

IS A -V A + (H -h)S c - S c x 


( 10 ) 


Let us introduce the following notations: b l = lS A - V\ and b 2 = 1S A - V 2 + (H - h)S c . 

Using these notations, the expressions for the pressure on the piston wall (5), (9), and (10) may be rewrit¬ 
ten as 


m ] RT l m 2 RT 2 m 2 RT 2 m 4 RT l 

P 1 b l + xS c ’ P 2 Z? 2 -x 5 c ’ b l + xS c ' P 4 b 2 — xS c ' ^ 

The resulting force acting on the rod with pistons is 

F rc = (Pl~ Pl + Pi~ Pi)S c . (12) 

The dynamic equation for the rod with pistons has the form 

m tc x = F rc (x) +F tc (x), (13) 

where F rs (x) is the velocity-dependent resistance force and m p = 2m ps + m rd . 
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Dynamic equation (13) should be complemented with the dynamic equations describing transport of the 
working medium’s mass through regenerators 

dm x dm 4 

—— = r(p l - p 2 ), — = —r(p A — p 3 ) (14) 

and the mass-balance equation 


m | + m 2 + m d 


m 0 1 , m 3 + m 4 + m d = m 02 . 


(15) 


Equations (13)—(15) form a closed system of nonlinear differential equations of the fourth order. It can 
be solved only by numerical methods (owing to nonlinearity), as is the case with the more simple models 
of external-heating engines [4, 5]. 

However, for such a model, it is reasonable to find an approximate solution for the dynamics that would 
enable analysis of a steady-state periodic process in a closed form. To develop such approximate model, we 
make the following assumptions: function F p (x) is linearized with regard to the central position of the rod’s 
displacement, i.e., the point x c = (H - /;)/2; continuously changing gas masses in hydraulic accumulators m ,, 
m 2 , m 3 , and m 4 are approximated with piecewise constant functions of the direct and inverse motion of a rod; 
and self-vibrations are assumed to be symmetric. 

Owing to the second and third assumptions, the expression for pressure values (11) for the direct motion 
of a rod may be presented as 


p iO) = 


b x + xS c 


Pi(x) = 


b 2 - xS, 


, Pi\(x) = 


-, p 4 (x) = 


b x + xS , 

where a x = m 1 RT 1 , a 2 = m 2 RT 2 , a 3 - m x RT 2 , and a 4 — m 1 RT l . 


b 2 -xS c ’ 


In accordance with (12), the resulting force may be written in this notation as 



a i 

b x + xS c 


a-, a , 

--— +--— 

b 2 -xS c b x +xS c 


_V 

b 2 -xS c J c ' 


(16) 


When the rod is located in the central position x c = (H - h)/2, (16) has the form 


pn 


H-h 
2 

, (H-h 


{a x + a 3 )/ 


pn 


' ^H-h„ 

b\ + —x—S c 


H-h 


- ( a 2 + a 4 )/ 


■ H-h „ 
b 2 - x~S c 


S c , 


(a l + n 3 )/ b ! + — —S c +(a 2 + a 4 )l\ b 2 


H-h. 


2 -, 


Sr- 


(17) 


Linearization of the function F p (x) taking into account (17) yields 

^pO) = F p (x c ) + (x-x c )F;(x c ). (18) 

We assume that the resistance force F rs behaves as the linear viscous-friction force F rs = kx. Substituting 
linearized expression (18) instead of the F p (x) force into (13), we arrive at a differential equation describing 
the rod’s motion from left to right, m p x + kx - (x - x c )F p (x c ) = F p (x c ). 

The equation describing the rod’s motion from the right to the left has a similar form. However, one may 
make an approximation that the term F p (x c ), while having the same absolute value, changes its sign to the 
opposite one. 
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Taking into account that F ' (x c ) is negative in 
both cases and introducing displacement with 
respect to the rod’s central position x c , x t = x - x c , the 
equation for the rod’s motion in both directions may 
be presented as 

m p x + kx l + |Fp(x c )|x 1 = T’ p (x c )sgnx 1 . (19) 

Equation (19) is identical to the known equation 
for self-vibrations of a system with negative dry fric¬ 
tion [8], so that its rigorous analytical solution may 
be found by the fitting method. It is reasonable to 
start with its transformation to the form 

Xj + 2kx x + COqXi = / p (x c ) sgnxj, (20) 

2 i 

where 2k = K/m p , co 0 = | F p (x c )|/m p and f p = F p (x c )/m p . 

A solution to Eq. (20) in a steady-state periodic 
regime has the following form: for the motion from 
the right to the left (0 <t< 7 / 2 ), 

Xj(?) = (x* + ^-\\e k '{coscot + —sincofl-^-|;(21) 
^ GV ^ W ' COo 



Xi(t) 




and for the motion from the left to the right (7/2 <t<T), 



r c t\] 

k . 

r c t\] 

COS 

CO \t- - 

-\—sin 



L l 2J] 

CO 

L l 2)\ 



( 22 ) 


2 , 2 


where to = ^ot 0 -k is the self-vibration frequency, 7 = 27t/co is the period of self-vibrations, 


xf = (/ p /co 0 z ) 


nk 


l+e m \/ 1-e 


nk 


is the amplitude value of self-vibrations corresponding to the right- 


nk 

(0 


nk 

to 


most position of a rod, and x^T/2) = -x* e - (/ p /co 0 ) |^1 + e j is the leftmost position of a rod (before 

the start of its motion to the right). 

Figure 2 shows a plot of the rod’s leftward motion in the steady-state regime x( t) = x,(?) + x c calculated 
according to (21) and (22) for the following values of parameters: H = 0.5 m, h = 0.03 m, / = 0.2 m, S A = 
0.031 m 2 , S c = 0.126 m 2 , p = 1.29 kg/m 3 is the density of the working medium (air), R p = 8.31 J/(mol kg), 
M = 0.029 kg/mol, m wm = 0.021 kg is the total mass of the working medium, m d = 0.0084 kg, 7, = 750 K is 
the heating temperature, 7 2 = 295 K is the cooling temperature, V t = V 4 = 0.0082 m 3 is the first and the fourth 
total volume of the liquid, V 2 = V 3 = 0.126 m 3 is the second and the third total volume of the liquid, m, d = 
10 kg, and m ps = 5 kg. 


Figure 3 shows the dependences of the vibration frequency to = ,J\F' p (x c )\/m p - (k /4m p ) on the heating 
temperature T l plotted for different values of the resulting mass m p = 20, 24, and 31 kg. 
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